Saturation in Deep Inelastic Scattering 
from AdS/CFT 



Lorenzo Cornalba a , Miguel S. Costa 6 

a Centro Studi e Ricerche E. Fermi 
Compendio Viminale, 1-00184, Roma 
Universita di Milano-Bicocca and INFN, sezione di Milano-Bicocca 
Piazza della Scienza 3, 1-20126 Milano, Italy 

b Departamento de Fisica e Centro de Fisica do Porto, 

Faculdade de Ciencias da Universidade do Porto, 
Rua do Campo Alegre 687, 4169-007 Porto, Portugal 

Lorenzo . CornalbaOmib . inf n. it , miguelcOf c .up .pt 



Abstract 

We analyze deep inelastic scattering at small Bjorken x, using the ap- 
proximate conformal invariance of QCD at high energies. Hard pomeron 
exchanges are resummed eikonally, restoring unitarity at large values of 
the phase shift in the dual AdS geometry. At weak coupling this phase is 
imaginary, corresponding to a black disk in AdS. In this saturated regime, 
cross sections exhibit geometric scaling and have a simple universal form, 
which we test against available experimental data for the proton structure 
function F 2 (x, Q 2 ). We predict, in particular, the dependence of the cross 
section on the scaling variable (Q/Q B ) 2 in the deeply saturated region, 
where Q s is the usual saturation scale. We find agreement with current 
data on F 2 in the kinematical region 0.5 < Q 2 < 10 GeV 2 , x < 10 -2 , 
with an average 6% accuracy. We conclude by discussing the relation of 
our approach with the commonly used dipole formalism. 

1 Introduction 

The high energy behavior of QCD is greatly simplified by the asymptotic weak- 
ness of the coupling and the approximate conformal invariance of the theory. Of 
great interest, in this respect, is the study of interaction processes in the Regge 
limit of high center of mass energy, with the other kinematical invariants kept 
fixed. This kinematical regime is, for instance, relevant to the analysis of deep 
inelastic scattering (DIS) experiments at fixed photon virtuality Q 2 in the limit 
of vanishing Bjorken x. 



Figure 1: Impact points p, x and p, x in H 3 , separated by a geodesic distance 
B. Dropping the radial coordinates p and p, we obtain the impact points x and 
x in the plane transverse to the high energy process R 2 , with impact parameter 
b = x — x. 



To the extent that QCD can be approximated by a conformal field theory, 
we must in general analyze CFT correlators of the form 



in the limit of large s — — (qi + q%) , and at fixed virtualities Qf = qf and 
momentum transfer t = — (51 + 473 ) 2 . In the high energy limit, the correlator 
([TJ is best analyzed in impact parameter space. The correct representation is 
suggested by the AdS/CFT duality [1], although let us stress that all of the 
results in this paper are purely based on simple implications of conformal sym- 
metry and could be derived also in the field theory language. Considering (JlJ 
as a high energy process in AdSs, the relevant transverse space is then the 
three dimensional hyperbolic space H3, holographically dual to the usual two- 
dimensional plane transverse to the high energy process described by (JT|), as 
shown in Figure [l] Representing four dimensionl vectors as x), with 

x ± lightcone variables and x a two dimensional transverse vector, we can pa- 
rameterize H3 using Poincare coordinates /?, x with metric p~ 2 [dp 2 + dx 2 ) and 
volume form p~ 3 dp d 2 x, with p > the distance to the holographic boundary 

of H 3 . Following the results in gJElll[51|5J[71[Elli[TU], we may write the 
impact parameter representation for the correlator ([TJ . Choosing, for simplicity 
of exposition, external scalar operators, it is given by 



where q is the transverse momentum transfer with — t = q 2 and where b is the 



(OM 2 (q 2 )0*(q 3 )0*(q A )) 



(1) 




(2) 
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usual impact parameter. The phase shift <5(s,b) is itself given by 

e 2i 5(s ,b) = J m m J | m m e2lA (S,B) (3) 

with A(S, B) the phase shift in AdS, which depends on the AdS energy squared 
and impact parameter S and B, according top] 

S = pps , 

cosh B= l±t±^. (4) 
2pp K > 

In particular, B is the geodesic distance between the points p, x and p, x in H 3 , 
with b = x — x. These represent the impact points of the operators 0\ and O2 
in the transverse space. Finally, the functions /$ are the radial wave functions 
for the scattering states. For scalar operators Oi, O2 of dimension Ai, A2 they 
are given by /, <x Q l p 2 K Al -2(Qip) for i = 1,3, and by f, oc Qip 2 K& 2 ^ 2 {.Qip) 
for i — 2,4 [TT]. We normalize the wavefunctions so that 

%h{ P )h(p)= I %h{p)k{P) = i- (5) 
P A J 

As shown in [5], the impact parameter representation ^ and (|3) approx- 
imates the conformal partial wave decomposition of the correlator (11) in the 
channel O1O2 — > 0\0* 2 , with intermediate states of conformal dimension and 
spin respectively given by y/S cosh(B /2) and v / 5'sinh(i?/2). In analogy with 
the usual results for scattering in flat space, we then expect that AdS unitarity 
implies [9J QH] 

ImA(5,B)>0, 

even though the phase shift 5(s, b) does not satisfy a simple unitarity constraint. 

We shall focus, for concreteness, on the very relevant and simple case of 
vanishing momentum transfer q = and equal virtualities for the incoming and 
outgoing states Q — Q\ = Q3 and Q = Q2 = Qi- It is then natural to construct, 
from the correlator Q, the following effective cross section 

S(s, Q,Q)= 2 J d 2 b Re (l - e M ^ b )) . 

Using ([3| , the cross section S can be conveniently written as 

^ hip) hip) J %hip)hiP) <?is,p,p) , (6) 

where we have defined the unintegrated cross sections 

o-(S) P> P) = J rf2b a ( s > P> P> b ) 1 ( 7 ) 
a( S ,p,p,b) = Re(l- e 2lA ^ B ) 



x We take the AdS quantities S and B to be dimensionless, measured in units of the AdS 
radius. 
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In this language, a(s, p, p, b) is the natural object which automatically satisfies 
the unitarity bound < a < 2 due to AdS unitarity ImA > 0. Moreover, for a 
black disk region we have a — ► 1, corresponding to a phase shift A with large 
imaginary part. 

In general, we cannot evaluate the integral over the impact parameter b. 
We may, on the other hand, use the relation Q between b and the AdS impact 
parameter B to rewrite the cross section a(s, p, p) in Q as 

cr(s, p, p) = 2npp f dB sinh B Re (l - e 2iA ^ p P' B A . (8) 

J\Mp/p)\ V ' 

It is now apparent that we are probing the phase A(S, B) for fixed S — spp and 
for B > |ln(p/p)|. Finally, note that the unintegrated cross sections satisfy, due 
to conformal invariance, non trivial relations under the transformation p — > p 2 / ' p 
with s — > s (yd 2 //? 2 ), b — > b (p/p), which leave invariant S and B. More precisely, 
a(s, p, p, b) is invariant whereas 

p 2 ( P 2 f _\ , 
P \ P P J 

The phase shift A (S, B) depends both on the number of colors N and on 
the 't Hooft coupling a s — a s N/ir of the theory. For large energy squared 
and impact parameter S and J5, the phase A will be dominated by the leading 
Regge pole of the planar diagrams of the theory [SJ [TD] , and will have a general 
representation of the forrrj^] 

A(S, B) = ^Jdv p{v) S^~ x fl iv (B) , (9) 

where the Regge spin j(y) and residue f3(v) depend implicitly only on the 't 
Hooft coupling a s and are even functions of v. The function fli U (B) computes 
radial Fourier transforms in H3, satisfies (Oh 3 + v 2 + l) £li U = and is given 
explicitly by 

o fn\- 1 vsinuB 

Whenever the AdS phase satisfies |A| -c 1, the full cross section is well 
approximated by a single Reggeon exchange, and we may write 

a(s,p,p,b) ~ 2ImA(5,.B). (10) 

In this case, the integral over the impact parameter b can be explicitly per- 
formed. In fact, using the Regge representation (Tol) for the phase shift, together 
witl£] 

1 -J? 

J :) / 



J d 2 b n iv (B) = t^pp 



2 In this paper, in order to have a uniform notation, we use slightly different conventions 
then in |9l I10|. In particular, A = — Trr t h ere and /3h er e = — '""Ahcre- 

3 In order to correctly compute the normalization of this Fourier transform, as well as of 
the ones in the sequel of the paper, it is safest to compute at non-zero momentum transfer 
and take the limit q — * 0. 
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coming from the integral representation |10j 

\ l+iv 



2 + / \p 2 + z 2 



s p 2 + (b-z) 
we may evaluate the cross section <r(s, p, p) to be 



1 — iv 

P 



a(s,p,p) ~ Im J dv ${v) (sppy^- 1 IH . (11) 

2 The Cross Section Deep into Saturation 

At fixed AdS energy squared S, the phase A(S, B) will in general vanish in the 
limit B — > oo. On the other hand, as we approach smaller and smaller impact 
parameters, A will in general grow and reach saturation at B ~ Bs(S'), where 
A is of order one. 

We will be mostly interested in weakly coupled gauge theories, where the 
phase is predominantely imaginary [12 In this case, saturation is reached aiQ 

2ImA(S',B s (S*)) ~ 1 . 

A typical plot of the saturation line in the (B, In S^J plane is given in Figure [2] 
In particular, for large S we have the linear relation 

B 8 (5) ~lu lnS+--- , (12) 

where • • • represents subleading terms in S. This can be shown, as customary 
[T7] . by approximating the integral in (JsJ) at the saddle point iB = j'(v)lnS. 
Saturation is then reached when the phase in ([9| vanishes at the saddle - i.e. 
when (1 + iv s ) B s = (j(^ s ) — 1) lnS". These conditions imply that 

u = -i j'(v a ) 

where the saturation saddle point v a is defined in terms of the Regge trajectory 
by 

(1 + iv s ) uj = j{v s ) - 1 . 
The cross section a(s,p,p) near saturation |ln(p/p)| > B s (spp) exhibits ge- 



ometric scaling [18 . More precisely, the integral (111 has a leading behavior 
given by 

a(s lPl p)~p 2 r -(i+-=)^ , (13) 



4 On the other hand, at strong coupling and for large impact parameters the phase shift is 
predominantly real and is given by the gravi-reggeon exchange in AdS. Studies of DIS in this 
regime include |14l . Saturation effects at strong coupling have also been analyzed in 115) . and 
a conjectured relation to black hole formation was put forward in |16| . 

5 Note that it is usually believed that, when ImA ~ 1, non-linear BK corrections to A 
of order N~ 4 due to fan diagrams also become relevant 1241 . As long as those corrections 
are negligible for impact parameters larger then the saturation line, they are irrelevant in the 
discussion which follows, since they will predominantly affect the phase shift in the black disk 
region. 



5 




Figure 2: Saturation line B S (S') in the B-ln(S) plane. As we increase S, the 
saturation line starts at a minimal value of S of the order of aj 1 and reaches the 
asymptotic linear behavior, shown with a dashed line, for large S. We extended 
the graph to the left of the B = axis symmetrically, drawing the mirror image 
of the saturation line. This is convenient since B = B(p,p,b) is invariant under 
p «-> p and we wish to show separately regions with p > p and p < p. 



where we have defined the scaling variable 

t 
P 2 



P I 2\ 

la \ S P ) 



(14) 



On the other hand, we are interested in the analysis of the cross section er(s, p, p) 
deep inside saturation, that is for 



\ln(p/p)\<B s (spp) 



(15) 



In this case, the integral ^ is dominated by the region B < B s , where we may 
replace a(s, p, p, b) ~ 1. This situation corresponds to a simple black disk in 
AdS, even though this is less transparent from the four dimensional perspective. 
We then obtain the approximate expression for the cross section 



°{s,P,p) - 2 Kpp 



npp 



dB smhB 

Mp/p)\ 

2 coshB s (s pp) - P - - 
P P 



(16) 



Moreover, when S — spp is large and we are in the linear regime (12 1, we have 
the simpler approximate expression 



<?(s,P,p) - npp 



{sppf + (spp) 



-u> _ P _ P 

p p 



(17) 
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where we neglect any subleading term in ( 12 1. Note that, for B s 3> |hi(p/ p) \ 1 
equation ( 17 1 is dominated by the first term and we obtain 



a(s,p,p) 



P 2 r 



to be contrasted with (13), valid near saturation. Finally note that, even in the 



deeply saturated region, the cross section cr(s,p, p) grows with s with a power 
law, violating the Froissart bound. Recall though that (17) has been derived 



by assuming an exact conformal symmetry and that, in a conformal theory, the 
Froissart bound is not relevant since there is no mass scale. 

For the sake of clarity, let us discuss a specific example, which is at the same 
time simple and instructive, since it contains most of the relevant physics. We 
will work with the maximally superconformal version of QCD, TV = 4 SYM with 
SU (N) gauge group, and we will consider the scalar operators 0\ — Tr (Z 2 ) 
and 2 = Tr (W 2 ) of dimension Aj = A 2 = 2, with Z and W two of the three 
complex adjoint scalars of the theory. To leading order in a s , we have the well 
known BFKL result 12 



1 



1 



and HHl 



„ , An tcLllll ~FT~ 

p{y) ~ i \6n A o? s 



v{l + v 2 y 

At vanishing In S, the integral ^ can be explicitly computed to be 

„-B 



(18) 



A(S=l,B) = ^a 2 s 



(6B 2 + 12B - Tr 2 ) 



sinh B 



-121n(l - 



6 

tanh B 



Li 2 (e- 2B ) 



In particular, we see that at B = we have 2ImA(S' = 1,B = 0) ~ 6.6 a 2 , 
which for a typical value of a s is well below saturation. At B = the saturation 
line starts for aslnS* ~ 1, as can be seen from the integral expression (|9| for 
the phase. The asymptotic linear regime (12) is reached for In 5 > 2/a s , with 
iv a ~ 0.26 and 

ui ~ 2.44 ct s . 

The above are clearly leading order results. However, as we shall explain in 
more detail in the next section, the experimental value of lo in DIS experiments 
is lower. For example, in the analysis of |19j one finds u> ~ 0.14, since the 
scaling variable t has the form (14) with 2uj/ (1 — lu) ~ 0.32. Therefore, as 
is well known, next to leading order corrections to the leading BFKL results 
(which also distinguish between QCD and its supersymmetric extensions) are 
important to match to experiment. 
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3 Deep Inelastic Scattering in QCD at Small x 



We now explore the phenomenological consequences of our results on deeply 
saturated cross sections for DIS in QCD at small Bjorken x. Throughout the 
discussion, we shall assume that we are working in the conformal setting, thus 
neglecting the running of the coupling constant and all quark masses. We will 
associate the scalar operators 0\ and O2 respectively to the photon and pro- 
ton. Note that, deep into saturation, the spin of the external particles plays 
a minor role, since amplitudes are dominated by the black disk region with 
cr(s,p, p, b) ~ 1. As usual, Q 2 is the photon virtuality and s ~ Q 2 /x. More- 
over, the scale Q will now represent a phenomenological parameter, related 
to the proton wavefunction, of the order of the relevant proton scales. The 
wave functions f\, fs and f± are localized respectively around p ~ Q^ 1 and 
Q^ 1 . Therefore, the total cross section E(s, Q,Q) in ^ can be approxi- 



mately computed using the saturated cross section a(s,p,p) in (16 1 whenever 

|ln (Q/Q) I < B s (s/QQ) . (19) 
Moreover, for large s/QQ, the saturation line B S (S') is approximately linear 



and a(s, p, p) is given by the simple expression ( 17 1. In this case, we may easily 
compute the radial integrals in ^ since, on purely dimensional grounds, we 
must have that 



dp 



hip) Hp) p c = Q~ c 7(0 



for some constants 7(C) of order unity, and similarely for the proton wave func- 
tions. Hence, in the deeply saturated regime at high s/QQ, we expect a rather 
simple form for the total cross section T,(s,Q,Q). Recalling that the cross 
section E is proportional to Q~ 2 F 2 , where F 2 {x,Q 2 ^ is the usual DIS proton 
structure function, we obtain that 



F 2 (x,( 



9. 

A 



Q_ 

xA 



Q_ 
xA 



- Q 


Q 


A 


— c — 
A 


A 





(20) 



where the constants c, c and the scales A, A are the only remenants of our lack 
of precise knowledge of the scattering radial wave functions. In particular, A, A 
will be of the same order as Q. Had we included the spin of the particles 
in the discussion, the parameters c, c, A, A would carry also this kincmatical 
information. The exponent lu is, on the other hand, universal and depends 
uniquely on the spin of the pomeron. Note that, since w < 1, the constants of 



order unity coming from the first two terms of (17 1 are essentially identical, and 
we may safely take 

A~Q 

in first ap pro ximation. 

As in (|17|, when B s (s/QQ) ^> |ln(Q/Q)| ^> 1, the cross section E is domi- 
nated by the first term in ( 20 1 and exhibits the geometric scaling 




Figure 3: Shown are the available measurements of F 2 (x,Q 2 ) in .21J, in the 
log 10 (Q)-log 10 (<3/x) plane, with energies measured in GeV. All points lie above 
the line x — 1. Shown is also the shaded region of points considered when 
analyzing ( p0| . It corresponds to the region delimited by the vertical line, setting 
Q > Qmin, the horizontal one, setting log w (Q/xQ) > 77, and the asymptotic 
linear saturation line. This is shown with a thicker line and is obtained by 
offsetting the graph in figure ^lby log 10 (Q) along the x = 1 line. 



where we define the scaling variable t as usual as 



91 

Qr 



Q 's = Q 



Recall that the power of 1 jx in the saturation scale Ql is observed experimen- 
tally, following [19';, to be given by 2w/(l-w) = 0.321 ± 0.056, so that 



0.138 ±0.021 



(22) 



Deeply into saturation, we predict that S evolves with r with the specific expo- 



nent in (21 1, which is uniquely fixed by the measurement of u> at the saturation 
scale Qf. This has to be contrasted with the behavior of £ near saturation 



following from (13 1, where the exponent of r is not fixed uniquely by u>. 
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We wish to test this prediction against the available experimental data on 
F2(x,Q 2 ). Measurements have been performed at values of x and Q 2 shown 
in Figure [3] as discussed in [20], which collects all available data from [2Tj . 
In this Figure, we present the data in the In Q-ln Q/x plane, with all energy 
scales measured in GeV from now on. These are the natural variables to discuss 



saturation, since they enter directly into ( 19 1. We will fit the available F2 values 



using ( 20 1 in its plausible region of validity. First of all, we will take 

Q ^* Qmin 

with Qmin ~ 1 GeV so that, in first approximation, the running of the cou- 
pling can be neglected. Secondly, we wish to choose points inside the saturation 



line (19 1. The exact determination of this line depends crucially on the phe- 
nomcnological parameter Q, and in turn on the strongly coupled dynamics of 
the proton. We expect the value of Q to be in the range of available scales - i.e. 
the QCD scale and the proton mass. Assuming that radial wave functions are 
localized around p ~ Q^ 1 and p ~ Q _1 , the saturation line in the lnQ-lnQ/x 
plane is then given by the saturation line for the AdS phase A shown in Figure 
[2J where we replace B and InS respectively by |lnQ/Q| and lnQ/(xQ). In 
practice, this ammounts to drawing the saturation line of Figure [2] onto Figure 
[3J offsetting the origin along the line x — 1 by InQ. We shall then take points 
with 

w In > In % , (Q ~ 0.2 - 1 GeV) . 

xQ Q 

Thirdly, we wish to consider data points with high values of Q/(Qx), so as to 
be into the linear regime of the saturation line. As explained in the previous 
section, the leading order BFKL analysis suggests that the linear regime starts 
around aiS > 2/a s , so that we shall take data with 

^ > 10V fa > 3) • 
Qx 

To proceed, let us choose Q m i n = 0.7 GeV, Q — 0.6 GeV and 77 = 3. We 
will show later that the main results are insensitive to this specific choice. The 
selected data is shown in the shaded region of Figure [3] We shall test our 
theoretical prediction against the real data in [21] as well as against the very 
accurate neural network interpolation to world DIS data in [22 . In particular, 
we shall minimiz e th e average square deviation of the data from the predi cted 



theoretical form (20 1. Since the parameters c, c and c/A 2 enter linearly in (20 1, 
minimization reduces to a linear system parameterized by the single parameter 

As a function of u>, the parameters c, c and c/A 2 are easily determined both 



for real as well as for simulated data. As is clear from (20), the parameter 
u controls the growth of F2 as 1/x increases. More precisely, at fixed Q, the 
coefficients c and lo determine the slope as well as the convexity of the function 
F2 in the experimental region of interest, shaded in Figure [3] Unfortunately, the 
relevant kinematics is on the boundary of the currently accessible experimental 
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0.1 0.12 0.14 0.16 CU 



Figure 4: Error functions for real and simulated data, plotted as a function of lu. 
For comparison, both functions have been normalized to 1 at their respective 
minima. The error for the simulated data exhibits a sharp minimum, while 
the corresponding error for real data is more insensitive to to, due to a large 
experimental uncertainty. We therefore choose to fix lu using the simulated 
data. The yellow stripe corresponds to the range u> = 0.1 38 ±0.021 obtained 
from geometric scaling. 

settings, resulting in data of relatively poor quality with large experimental 
uncertainty. This is reflected in the fact that the error function, although it 
presents a minimum for lu ~ 0.136, is essentially constant in the range 0.1 to 
0.17 plotted in Figure [4] On the other hand, as shown in the same figure, if we 
use the more accurate simulated function F 2 computed at the same values of x 
and Q 2 available in the real data [22j . we obtain a rather sharp minimum for 
the error function at 



Therefore, from now on, we shall determine the optimal value of lu using the 
simulated data only. At this point we wish to emphasize that this value of 
lu, obtained from data inside the saturated kinematical region, is within the 
experimental range lu = 0.138 ± 0.021, obtained independently from geometric 
scaling. The values for the other relevant parameters can be determined to be 
A ~ 1.0 GeV, c ~ 0.13, c ~ 0.14 for the real data and A ~ 1.0 GeV, c ~ 0.11, 
c ~ 0.08 for the simulated one. 

The real data is presented in Figures [5] and [6] where we show the data 



together with the theoretical curves from ( |20| . We plot F 2 /Q as a function of 
log 10 Q, and each graph contains data points with values of log 10 (<5/x) in the 
range Y ±0.1, for 2 < Y < 5 in increments of 0.2. Theoretical curves are shown 
in red both for the minimal and maximal values of log 10 (Q/izi). Finally, the 
shaded area corresponds to the region delimited by the choice of parameters 
Q,Vi Qmin, as also shown in Figure [3] Analogously, Figures [7] and [8] show the 



lu ~ 0.126 . 
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Figure 5: Real data F2/Q as a function of Q. Each graph contains data points 
with values of log 10 (Q/x) in the range Y ± 0.1, for 3.6 < Y < 5 in increments 
of 0.2. Theoretical curves are shown in red both for the minimal and maximal 
values of log 10 (Q/a;). The shaded area corresponds to the region delimited by 
the choice of parameters Q m i n = 0.7 GeV, Q = 0.6 GeV and 77 = 3, as also 
shown in Figure [3] 
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0.5 1.0 2.0 4.0 0.5 1.0 2.0 4.0 



Q [GeV] 

Figure 6: Real data F2/Q as a function of Q. Each graph contains data points 
with values of log 10 (Q/x) in the range Y ± 0.1, for 2 < Y < 3.4 in increments 
of 0.2. Theoretical curves are shown in red both for the minimal and maximal 
values of log 10 (Q/a;). The shaded area corresponds to the region delimited by 
the choice of parameters Q m i n = 0.7 GeV, Q = 0.6 GeV and rj = 3, as also 
shown in Figure [3] 
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Figure 7: Same as Figure [5] for the simulated function F% computed at evenly 
spaced values of Q for fixed Y = log 10 (Q/x). 
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Figure 8: Same as Figure [6] for the simulated function F2 computed at evenly 
spaced values of Q for fixed Y — log 10 (Q/x). 
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simulated data. 

Let us note that the AdS black disk form of the structure function given 
in (20 1, with the above choice of parameters, approximates the available real 
data with an average 6% accuracy in the rather large region of parameter space 
0.5 < Q 2 < 10 and x < 10~ 2 . 



Table 1: Number of experimentally available data points n and predicted value 
of lo for different values of Q and <2 mm - 



Q 


Qmin 


V 


n 


LO 


0.3 


0.7 


3 


58 


0.104 


0.3 


1 


3 


23 


0.090 


0.6 


0.7 


3 


138 


0.126 


0.6 


1 


3 


104 


0.130 


1 


0.7 


3 


200 


0.141 


1 


1 


3 


171 


0.152 



Due to uncertainty on the precise location of the saturation line, we have 
repeted the analysis with different values of Q, Q mm and 77, to test the robustness 
of the predicted value for lo. Within the range 0.7 < Q m [ n < 1 and 0.3 < Q < 1, 
the fitted value for lo varies from 0.090 to 0.152, as shown in Table [TJ thus 
mostly within the predicted range (the optimal value of lo is rather insensitive 
to the choiche of 77 > 3 which we keep fixed). Note that, although the first 
two entries of Table [l] are outside the predicted range, they are based on a very 
small number of data points. 

As already stressed, available data is on the boundary of the deeply saturated 
region, and one would need to reach higher energies in order to better test these 
predictions. Possibly, future data from LHC will be of use to confirm the above 
results. 



4 Relation to the Dipole Formalism 

We will conclude this paper by discussing the relation between the above results 
and the dipole formalism |23j . which is usually employed in the analysis of 
saturation effects. In this context, it is costumary to analyze the so-called 
dipole-dipole cross section <7dd(s, r, r, b) instead of a(s, p, p, b), with the full 
cross section S(s, Q, Q) given by integrals over the dipole transverse orientations 
r and r 

W(r) <7 DD ( S ,r,f) W(v) , (23) 
cr D D(s,r,r)= / d 2 h ct dd (s, r, r, b) , 
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Figure 9: Relation to dipole formalism. While the dipolc dipolc cross section 
cdd(s, r, r, y — y) depends on four points in R 2 , the cross section a(s, p, p, x— x), 
which is the one constrained by unitarity, depends only on two points in H3. 
The role of the dipole vectors r and r is now played by the radial coordinates p 
and p. 

where W(r) and W(f) are the so-called dipole impact factors. 

Let us first note that, although the dipole formalism is quite useful due to 
its intuitive physical description of the high energy process and of the linear 
BFKL and nonlinear BK evolutions [23], it is not well suited for the discussion 
of unitarization, since the natural object which satisfies the unitarity constraint 
< a < 2 is a(s,p,p, b), instead of <tdd(s, r, r, b). This fact is quite clear in 
gauge theories which are exactly conformal, like AT = 4 super Yang-Mills, where 
the dipole formalism can still be applied (as well as the BK equation, which is 
explicitly conformally invariant). In this case, the theory has no asymptotic 
states or an S-matrix to which to apply the usual unitarity constraints. More- 
over, even in a confining theory like QCD, which possesses asymptotic states, 
the dipole state is not a single particle state at infinity, and therefore does not 
enter in a usual S-matrix element. In fact, in the standard discussions of DIS 
at small x, the dipole picture is often used to describe the wave function of an 
off-shell spacelike photon. 

At any rate, in order to make contact with the standard litterature, we will 
briefly analyze, in what follows, the above expressions in the unsaturated regime 
of small AdS phase shift |A| < 1. 

Let us first analyze the impact factors W(r) and W(r), leaving to the sec- 
ond part of this section the discussion on 0dd(s, r, r, b) and on saturation in 
the context of the dipole formalism. We recall the BFKL representation of A 
analyzed in [TO]. More precisely, to leading order in the coupling we have that 

~^^d 2 yd 2 y VF(p,r,x-y) F(r,r,y-y) W{p, f,R-y), 
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with b = x— x. The impact factor W depends on the point p, x in H3 and on the 
two intermediate points y ± r/2 on the boundary of H3, as represented in Figure 
[9J Similar comments apply to the impact factor W. Moreover, F(r,r,y — y) is 
the leading order two-gluon exchange kernel from y±r/2 to y±r/2. Integrating 
against 



d 2 h 



%h{p)h{p) 



%f2{p)fi(p) 
P A 



obtain an expression of the form ( 23 1 , where 



and using the approximate relation p0| valid in the small phase regime, we 

(24) 



<t dd (s, r, r) ~ cc 2 Re / d 2 w F(r, r, w) , 

W(t)=- c J ^h{p)h{p) J d 2 w W(p, ri w) , 

and similarely for W(r). The constants c, c are fixed by the normalization 
conditions 

analogous to ([5]). 

Let us discuss the impact factor W in detail. As shown in [10], conformal 
invariance highly constrains W(p, r , w) to be a function of the unique cross ratio 



P' 

which can be conveniently written in the following integral representatiorj^] 
1 



64tt 5 



x / d 



dv v A (\ + u z 



r(i -iu) 



W(u)x 



(25) 



>-§) 2 (z+i) 2 , 



p 2 + (w - z) J 

where the transforms W(v) and W(v) determine the Regge residue (3{v) to be 

Mu) = — WM tanh(7r^/2) WM . 

Moving to momentum space in the transverse E 2 plane by integrating against 
J d 2 w we obtain 



1 

32tt 3 



dv {l + v 2 ) W(y) 



vc(i -c) 



C(l-C )r 2 
P 2 



2\ -2 



(26) 



6 In the notation of |10| . W(u) is given by V(v) /V m i n (u, 1), with V(v) the impact factor for 
the full amplitude and with V min (u, 1) = F ( 2Al ~ 1+ ^ ) T ( 2Al " 2 '~" / j /(r(Ai)r(Ai - 1)), 
where Ai is the dimension of the external operator 0\. Similarely for W. 
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where C is the Feynman parameter related to the denominators in the last 
parenthesis of (251. 

For concreteness, let us return to the specific example already discussed in 
section^ with (3(v) given by ( Jl8). We have that W(v) = 8ir 2 a s / (l + v 2 ). Then 

W(r) can be computed from (|26 L since the v integral fixes p = |r| (1 — C). 
After integrating against c" 1 Jdp P~ 3 fif3, we obtain 



W(r) 



o C(i-C) 



/i(yClW)|r|)/ 3 (yC(W)|i 



with c = a s /2. Using the fact that fi(p) = \/2Qip 2 K n (Qip) for i = 1,3, we 
have just obtained the usual expression in terms of dipole wave functions [2"3"| . 

To conclude this section, let us recall the form of the dipole-dipole cross 
section <7DD(s>r,r), in the unsaturated case (24). It is given by the usual BFKL 
kernel 



o-DD(s,r,r) 



2irN 2 



16-K 3 a 2 

dv t: ~ a 2 ( s l r ll r D 
(1 + ^) 



which should be compared with (111, with /3(z/) given by (18). Recall also 
that the above expression is the zero momentum contribution to the full BFKL 
expression for ctdd (s,r,r,b) given by [T3] 



cr DD (a, r, r, b) ~ ^N 2 \ dv v a ^ ( s l r ll r l) jW * ( r > r > b ) 



(27) 



where 



16^ r(i-^) r(i^) 5 



(i + ^) 2 r(i^) 2 r(i + »i/) ' 



and where % v (r,f, b) is the 2-dimensional conformal partial wave of spin 
and conformal dimension 1 +iv at the four points (b ± r)/2, (— b ± r)/2. Due 
to transverse conformal invariance, % v depends uniquely on the cross-ratio 
combinations 



zz 



(b- 



) 2 (b + |-|) 2 ' 



zz 



and is given explicitly by 

2i„ (r, f , b) = {~z) h {-z) h F {h, h, 2/i, z) F (ft, h, 2h, z) , 
h = — - — , 



where F is the hypergeometric function 2 -Pi. 
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The expression (27) should be confronted with a (s, p, p,h) derived from 
equations (10 1 and (|9j in the limit of small AdS phase shift |A| <C 1, 

0" (s, P, P, b) ~ / dv v Im („) S^- 1 . (28) 



2ir 2 N 2 J *~ x " ' " sinhB 

Given the similarity of the expressions for <tdd and a, we may try to follow 
the program of section 2 and consider the saturation region in the transverse 
impact parameter b-space where, at fixed energy s and dipole orientations r,r, 
the cross section odd (s, r, r, b) becomes greater then unity. On the other hand, 
this program cannot be carried out in general. In fact, while the cross section 
a (s, p, p, b) depends, aside from energy, on a single conformal cross ratio B, 
the cross section <tdd (s, r, r, b) depends on energy and on two cross ratios z, z 
which reflect the orientations of the two dipoles (the scalar quantities p, p are 
replaced by the dipole transverse vectors r,r). Moreover, both B and In S 
enter exponentially in ( 28 ) and allow for a simple saddle approximation of the 



integral and for the determination of the saturation line as in ( 12 1. On the other 



hand, the dependence of the integrand in (27 1 on z, z is now highly non-trivial, 
since it involves not only the norms but also the orientations of r, r, b. Most 
importantly, it does not allow for a simple approximation of the integral at a 
saddle point and a simple determination of the saturation line. 

The analysis of the saturation region for odd can be carried out only in the 
limit |r| , |r| <C |b|. In fact, in this limit, one may use the operator product 
expansion and obtain 

7^(r,r,b)c (29) 

of the simple exponential form. Similarly, i n th e same limit p,p <C |b| 
B ~ In (b 2 /pp) and one may substitute, in (28) 



has 



sinh B 



obtaining an expression analogous to ( 29 ) . From here on we may follow the usual 



steps reviewed in section 2 to determine the saturation radius, which is given 
in general by ln(b 2 /|r||r|) ~ wln(s|r||r|), leading to a black disk cross-section 
given by 7rb 2 ~ 7r|r||r|(s|r||r|) w . 

With this procedure we recover an expression analogous to the first term of 
(17) just as easily in the dipole formalism. On the other hand, as we already 



pointed out, the usual saddling argument works in (28 1 for generic values of 



p,p,h. This fact allows to determine the other terms in equation (17), which 



could not have been deduced in the dipole language. Note that the extra three 



terms in (17), and correspondingly in (20), are crucial in order to have a qual- 



itatively good fit for the relevant F2 data at hand. In fact, the expression for 
F2/Q exhibits a non-trivial dependence on Q at fixed Q/x, which is qualitatively 



correctly captured by the third and fourth terms in ( 20 1 proportional to 

_Q _ A 
A Q ' 
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These terms give a concave behavior with maximum at Q ~ A, which is a clear 
feature of the F 2 data, as can be seen from the plots in Figures [5] [6] [7] and [8] A 
pure term of the form (Q/x)^ , as could be determined by the above arguments 
also in the usual dipole formalism, is clearly insufficient to reproduce the Q 
dependence at fixed Q/x inside the saturation region. 

Acknowledgments 

We would like to thank Joao Viana Lopes and Giulia Galbiati for support 
with the data analysis, Andrea Banfi for discussions and Markus Diehl for 
carefully reading the manuscript. LC is funded by the Museo Storico della 
Fisica e Centro Studi e Ricerche "Enrico Fermi" and is partially funded by 
INFN, by the MIUR-PRIN contract 2005-024045-002 and by the EU contract 
MRTN-CT-2004-005104. MC is partially funded by the FCT-CERN grant 
POCI/FP/63904/2005. Centro de Fisica do Porto is partially funded by FCT 
through the POCI program. 

References 

[1] J. M. Maldacena, The large N limit of superconformal field the- 
ories and supergravity, Adv. Theor. Math. Phys. 2 (1998) 231, 
| JarXi"v:hep-th/9711200p . 

[2] J. Polchinski and M. J. Strassler, Hard Scattering and Gauge/String Dual- 
ity, Phys. Rev. Lett. 88 (2002) 031601, [ |arXiv:hep-th/0109174| ] . 

[3] R. C. Brower, J. Polchinski, M. J. Strassler and C. I. Tan, The Pomeron 
and Gauge/String Duality, [ |arXiv:h"e p-th/0603115|] . 

[4] L. Cornalba, M. S. Costa, J. Penedones and R. Schiappa, Eikonal approxi- 
mation in AdS/CFT: From shock waves to four-point functions, JHEP 08 
(2007) 019, | |arXiv : hep-th/061 1 122] | . 

[5] L. Cornalba, M. S. Costa, J. Penedones and R. Schiappa, Eikonal Approx- 
imation in AdS/CFT: Conformal Partial-Waves and Finite N Four-Point 
Functions, Nucl. Phys. B 767, 327 (2007) QarXiv : hep-th/061 1 123] ] . 

[6] L. Cornalba, M. S. Costa, J. Penedones, Eikonal Approximation in 
AdS/ CFT: Resumming the Gravitational Loop Expansion, JHEP 09 (2007) 
037, arXiv: 0707. 0120 [hep-th] . 

[7] R. C. Brower, M. J. Strassler and C. I. Tan, On the Eikonal Approximation 
in AdS Space, arXiv : 0707 . 2408 [hep-th] . 

[8] R. C. Brower, M. J. Strassler and C. I. Tan, On The Pomeron at Large 't 
Hooft Coupling, arXiv : 0710 . 4378 [hep-th]. 



21 



[9] L. Cornalba, Eikonal Methods in AdS/CFT: Regge Theory and Multi- 
Reggeon Exchange, arXiv: 0710. 5480 [hep-th] . 

[10] L. Cornalba, M. S. Costa, J. Penedones, Eikonal Methods in AdS/CFT: 
BFKL Pomeron at Weak Coupling, arXiv: 0801 . 3002 [hep-th]. 

[11] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Gauge theory cor- 
relators from non-critical string theory, Phys. Lett. B 428, 105 (1998) 
| ]arXiv:hep-th/9802109p . 

E. Witten, Anti-de Sitter space and holography, Adv. Theor. Math. Phys. 
2, 253 (1998) [ arXiv : hep-th/9802150"l l . 

[12] V. S. Fadin, E. A. Kuraev and L. N. Lipatov, On The Pomeranchuk Sin- 
gularity In Asymptotically Free Theories, Phys. Lett. B60 (1975) 50-52. 

E. A. Kuraev, L. N. Lipatov and V. S. Fadin, The Pomeranchuk Singularity 
in Nonabelian Gauge Theories, Sov. Phys. JETP 45 (1977) 199-204. 

Ya. Ya. Balitsky and L. N. Lipatov, The Pomeranchuk Singularity in Quan- 
tum Chromodynamics, Sov. J. Nucl. Phys. 28 (1978) 822-829. 

[13] L. N. Lipatov, Small-x physics in perturbative QCD, Phys. Rept. 286, 131 
(1997) QarXiv : hep-ph/96 10276] ] . 

[14] J. Polchinski and M. J. Strassler, Deep inelastic scattering and gauge/string 
duality, JHEP 0305 (2003) 012, [ ]arXiv : hep-th/02092li"] ] . 

C. A. Ballon Bayona, H. Boschi-Filho and N. R. F. Braga, Deep inelastic 
scattering from gauge string duality in the soft wall model, JHEP 0803 
(2008) 064, arXiv: 0711. 0221 [hep-th]. 

C. A. Ballon Bayona, H. Boschi-Filho and N. R. F. Braga, Deep inelas- 
tic structure functions from supergravity at small x, arXiv: 0712. 3530 
[hep-th] . 

[15] Y. Hatta, E. Iancu and A. H. Mueller, Deep inelastic scattering at strong 
coupling from gauge/string duality: the saturation line, JHEP 0801 (2008) 
026, arXiv: 0710. 2148 [hep-th]. 

[16] L. Alvarez-Gaume, C. Gomez and M. A. Vazquez-Mozo, Scaling 
Phenomena in Gravity from QCD, Phys. Lett. B 649 (2007) 478, 
[ ]arXiv:hep-th/0611312p . 

L. Alvarez-Gaume, C. Gomez, A. S. Vera, A. Tavanfar and M. A. Vazquez- 
Mozo, Critical gravitational collapse: towards a holographic understanding 
of the Regge region, arXiv : 0804 . 1464 [hep-th] . 

[17] L. V. Gribov, E. M. Levin and M. G. Ryskin, Semihard Processes In QCD, 
Phys. Rept. 100 (1983) 1. 

E. Iancu and L. D. McLerran, Saturation and universality in QCD at small 
x, Phys. Lett. B510 (2001) 145 [ ]arXiv : hep-ph/010303"2] l . 



22 



A. H. Mueller, Parton saturation at small x and in large nuclei, Nucl. Phys. 
B558 (1999) 285 QarXiv : hep-ph/99044d4| ] . 

A. H. Mueller and D. N. Triantafyllopoulos, The energy depen- 
dence of the saturation momentum, Nucl. Phys. B640 (2002) 331 
| ]arXiv:hep-ph/0205167l ] . 

S. Munier and R. Peschanski, Geometric scaling as traveling waves, Phys. 
Rev. Lett. 91 (2003) 232001 [ ]arXiv : hep-ph/0309177l l . 

S. Munier and R. Peschanski, Traveling wave fronts and the transition to 
saturation, Phys. Rev. D69 (2004) 034008 [ |arXiv : hep-ph/0310357l ] . 

S. Munier and R. Peschanski, Universality and tree structure of high energy 
QCD, Phys. Rev. D70 (2004) 077503 QarXiv : hep-ph/0401215| | . 

[18] A.M. Stasto, K. Golec-Biernat and J. Kwiecinski, Geometric scaling for 
the total gamma* p cross-section in the low x region, Phys. Rev. Lett. 86 
(2001) 596 arXiv:hep-ph/0007192 

[19] F. Gelis, R. B. Peschanski, G. Soyez and L. Schoeffel, Systematics of geo- 
metric scaling, Phys. Lett. B 647, 376 (2007) | )arXiv : hep-ph/06 1 0435} ] . 

[20] L. Del Debbio, S. Forte, J. I. Latorre, A . Piccione and J. Rojo [NNPDF 
Collaboration], JHEP 0503, 080 (2005) |arXiv:hep-ph/0501067] . 

[21] M. Arneodo et al. [New Muon Collaboration], Nucl. Phys. B 483 (1997) 3. 

A. C. Benvcnuti et al. [BCDMS Collaboration], Phys. Lett. B 223 (1989) 
485. 

M. R. Adams et al. [E665 Collaboration], Phys. Rev. D 54 (1996) 3006. 
M. Derrick et al. [ZEUS Collaboration], Z. Phys. C 72 (1996) 399. 
J. Breitweg et al. [ZEUS Collaboration], Phys. Lett. B 407 (1997) 432. 
J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 7 (1999) 609. 
S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 21 (2001) 443. 
J. Breitweg et al. [ZEUS Collaboration], Phys. Lett. B 487 (2000) 53. 
C. Adloff et al. [HI Collaboration], Nucl. Phys. B 497 (1997) 3. 
C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 13 (2000) 609. 
C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 21 (2001) 33. 
C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 19 (2001) 269. 
C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 30 (2003) 1. 

[22] S. Forte, L. Garrido, J. I . Latorre and A. Piccione, JHEP 0205, 062 (2002) 
arXiv:hep-ph/0204232 . h ttp: / /sophia.ecm.ub.es/nnpdf/| 



23 



[23] A.H. Mueller, Soft Gluons In The Infinite Momentum Wave Function And 
The Bfkl Pomeron, Nucl. Phys. B 415 (1994) 373 

A.H. Mueller and B. Patcl, Single And Double Bfkl Pomeron Exchange 
And A Dipole Picture Of High-Energy Hard Processes, Nucl. Phys. B 425 
(1994) arXiv:hep-ph/9403256 

A.H. Mueller, Unitarity and the BFKL pomeron, Nucl. Phys. B 437 (1995) 
107 arXiv:hep-ph/9408245 

[24] A. H. Mueller, Parton saturation: An overview, arXiv:hep-ph/01 11244. 



24 



